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pared with the prediction of the parton model. Appreciable deviation is observed 
■ J . only at low values of the Bjorken scaling, which correspond to the extremely deep 

VI . inelastic scattering. 
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Abstract: The proton structure function is re-deduced from the data of deep in- 
elastic electron-proton scattering after enhanced correction that is made due to the 
multiple scattering effect. The Glauber approach is used to account for the multiple 
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1. Introduction 

A dynamical understanding of quark substructure had its origin from experiments 
on deep inelastic lepton-nucleon scattering [1]. These showed that the complicated 
production of many hadrons in such a collision could be simply interpreted as quasi- 
elastic scattering of the lepton by point like particles or quarks. For example the 
linear relation between the total inelastic cross section cr*°*(e — p)and the center of 
mass energy E may be explained as simple elastic scattering with point particle 
which depends on the coupling constant and the phase space factors. It is well 
known that at high momentum transfer, the elastic form factor is very small and 
the inelastic scattering of the incident electron is much more probable than the 
elastic scattering. In this case there must be an extra variable x (the Bjorken scaling 
variable) which relates the momentum transfer q to the energy transfer u by the 

2 

relation x = , < x < 1 . The case of x = 1 corresponds to the elastic 
scattering, and it tends to zero as the reaction goes toward the deep inelastic. The 
differential cross section [2] is then written as, 

-[F2(g^ z/) cos^ - + —Fr{q\ v) sin^ -] (1.1) 



dq^dv q^ ' ' 2 M ' ' T 

Where the factor represents the square of the scattering amplitude of the electron 
by a point charged particle (the parton or quark). F^ and F2 are the hadron structure 
functions used to correct for the composite structure of the target, i?and E' being the 
incident and scattered electron energies. The factor cos^ (|)stands for the square of 
the rotational matrix element for the electric non-spin flip part of the scattering of the 
electron by a charged point particle. The factor sm^(|) is the corresponding factor 
for the spin-flip magnetic interaction between the magnetic moments of the electron 
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and target particle. The problem arises now revolves about the factor 47rQ;^/g^ , 
which assumes a simple scattering of the electron by a Coulomb field of a charged 
point particle. And since the nucleon bag contains 3- valence quarks rather than the 
sea quark- antiquark pairs, so we expect multiple collision inside this bag instead of 
a single one. Also the presence of more than one quark in the system may make 
shadowing on the others during their collision. Of course these will affect the shape 
of the differential cross section and hence the structure functions Fi and F2. The 
present work aims to consider the effect of multiple scattering by using the Glauber 
[3] approach and then re-deduce the form of the structure function and compare 
the results with the previously published data at different regions of the Feynman 
variable x. 

2. The Multiple Scattering Approach 

Let us consider a bound system of three point particles (valence quarks) forming a 
proton as a spherical bag. The sea quark- antiquark pairs are considered as electric 
dipoles in their plasma phase. The problem is treated in the laboratory system in 
which its origin coincides with the center of the spherical bag. Assuming a collision 
axis in the z- direction, an incident electron with position vector and impact param- 
eter will scatter with a momentum transfer . The scattering amplitude is given by 

[4], 

fmsiq) = ^ / d^bexp{i q .b) {s,} |cxp(zS(fe, {s,}) - 1| ^{s,}) (2.1) 

Where, k' is the momentum of the scattered electron. In case of inelastic scat- 
tering, only a fraction x of the momentum of the incident electron will transfer to the 
quark system of the proton, leaving a momentum fraction {1 — x) to the scattered 
electron, is the total phase shift function due to all the constituent valence quarks 
and the associated dipoles of the system having position vectors {rj} with projections 
{si} on the impact parameter plane. The suffix i runs over all the flavors forming the 
proton bag. is the total wave function of the proton system. In Quantum Chromo- 
dynamics [5,6,7], the force between quarks is mediated by the exchange of massless 
spin 1, gluons, which are the carrier of the color force. When QCD is dimensionally 
reduces to 1 + 1 dimensions [8,9,10], the gauge fields may be eliminated using their 
equations of motion. What remains is the effectively long-range quark-quark force, 
which is given by a linear two body potential [11]. This linear potential can be un- 
derstood in other ways. Since gluons are massless, their propagator in momentum 
space is 1/q'^. This is the same as the photon propagator, which corresponds to the 
\r-r'\ Coulomb potential. An appropriate wave function for the valance quarks inside 
the proton should satisfy certain conditions. According to the boundary condition 
of the problem, the quarks are considered as free particles in the central region of 
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the nucleoli and they are tightly bounded at large distance. Consequently, the nu- 
cleon space is classified into two regions. The central part, where s < Si < sq , the 
quark-wave function is taken as plane wave. And in the peripheral region Sj > sq, the 
quark wave function is approximated to harmonic oscillator wave function. Here Sgis 
a parameter which specifies the boundary of the central region. In the approximation 
of the straight line trajectories [3], the total phase shift function may be expressed 
as the sum of the phase shift functions due to the 3- valence quarks and N- dipoles, 
as scattering centers, 

dip 

{bi}) ^T^Ub - bi) + j:Ub - bi) (2.2) 

i i 

is the phase shift due to the ith valence quark, and ^f^^is the phase shift 
function due to the i*'* electric dipole. Defining a profile function as r(6, {bi}), 

r(6,{6,}) = l-exp[iS(6,{6,})] (2.3) 
Which could be expressed in terms of the two-body profile function Fj, 

m {h}) = 1 - n[i - Ub - h)] n[i - (& - m 

3 N 3 

= (E + E rf ) - (E r^r,- + E ^f'^f + ■■■ (2-4) 

N 3 

+(-1)^-1 nrfnr. 

j i 

The terms of the first bracket in the R.H.S. of Eq.(2.4) represent the single 
scattering or the impulse approximation, the second bracket is the correction due to 
the binary collisions and finally the last term is the correction due to the collision 
with all the constituent of the proton bag, i.e. the coUision at which all the three 
valence quarks (uud) and the N- pairs contribute the reaction. Assuming that all 
the quark dipoles are similar and equally probable, then Eq. (2.4) becomes, 

r(6, {bi}) = (Er, + NV^'p) - (^r,r, + (r-^-^)^ (2.5) 

+ ---+(-i)^-i(r'^*p)^nr, 

i 

At sufficiently high energy, it is expected that the number N of created is much 
greater than the number of valence quarks. And if this number is linearly propor- 
tional to the center of mass energy Ecm of the reaction, then the terms of Eq. (2.6) 
are energy dependent, and read as. 
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r{b, {k}) E^rnT'-P - ^V '^")' + • • • (2.6) 

So that the first term in Eq. (2.6) is hnearly dependent on E„n while the second 
term that represents the double scattering is quadratically dependent. This term is 
not important even at high energy, since the electromagnetic fine structure constant 
a appears of second order through the profile function that depresses the term as a 
whole. The result of Eq.(2.6) is very important, since it explains clearly the linear 
proportionality of the e-p total cross-section with the center of mass energy. 

The problem now is reduced to the two-body collision, in other words, we have 
to formulate, only, for the simple two-body function and in terms of the scattering 
potential, 

^jib - bj) = 1 - ^^P[- 2A:(l-x) i-oo ~ ^^"^^^ ^^-^^ 
Assuming a Coulomb like potential field with the form, 

U — — exp(— ar) and f— {b — bj)+ ^ (2.8) 

r,(6 - h) = 1 - ^M^kir^^^oHb' + b] - 2bbj COS <f>j))] (2.9) 

Ko{b'^+b'j—2bb.j coscpj) is the modified Bessel function or the so called Mc-Donald 
function of order zero. Uo, being the strength of the field. For the two fiavors u and 
d of the valence quarks having fractional charges 2/3 and —1/2 respectively, we get 



Kib - bj) = 1 - e^P[-^^(Yr^^°(^«(^' + ^? - 2&&, cos0,))] (2.10) 
And 



^dib - bj) = 1 - exp[- 



2fc(l 



X 



-Ko{--a{b^ + b] - 2bbj cos(/),))] 



(2.11) 



To calculate the profile function Tf^(b — 6j)due to the scattering from a single 
dipole, let us consider that a distance £ separates the two quarks of the dipole. The 
potential field due to the dipole is given by. 



U, 



A 



+ 



A 



r-raip + I/2 f - r^ip -1/2 



(2.12) 



And putting, r—b + z, Tdip—bdip + z , i—bt -\- z where z is the projection along 
the collision axis and b is the projection in the impact parameter plane. So that the 
phase shift function due to a dipole is. 



4 



JHEP 2.0 DRAFT (draft, published, notoc) 



JHEP, February 1, 2008, 12:49 



and the total phase shift due to the N-dipoles is, 
Then Eq.(2.1) becomes, 



fmsiQ,^) = / d' b exp(i q ■ h) |r(6, {s,}\ (2.15) 

3. Results and Discussion 

The model predictions are tested with the experimental data of SLAC-E-049, SLAC- 
E-061, and SLAC-E-089 [12] that concerns the deep inelastic scattering in the lab 
momentum ranges {3.748 - 24.503}. A spherical proton bag is assumed with root 
mean square radius Rrms = 1-2 fm. Since the quarks are almost considered inde- 
pendent, then the proton spatial wave function is taken as, "i/p = Y[u,u,d^i ) with $j 
as a plane wave for s < So and $i is a harmonic oscillator wave function for s > So- 
The parameter So is adjusted to fit with the e — p inelastic cross-section at the cor- 
responding center of mass energy. It is found that Sq is ^Rrms at the considered 
energy. The total profile function of the e-p scattering is calculated with a weight 
factor corresponding to the quark probability density inside the proton bag. 

Fig.(l) shows the individual terms of the average profile function corresponding 
to the single, double and the triple scattering. Of course the first term is the most 
important one allover the impact parameter range. The double scattering term 
represents about 0.02 from the first one at the central region. This ratio is greater 
than the reciprocal of the fine structure constant The central region is populated 
with high density of states, so there will be a good chance for the double collision 
to occur. The values of the triple scattering may be neglected. In other words, the 
shadowing has significant effect only up to the double scattering. The resultant profile 
function is plotted in Fig. (2), which shows slow decrease with the impact parameter b, 
reficcting the Coulomb behavior of the assumed scattering potential. The ratio of the 
imaginary to the real part of the profile function is an important coefficient. It is just 
a measure of the degree of inelasticity of the reaction. As demonstrated in Fig. (3), 
the imaginary to the real part of the coefficient shows that the central region is 
characterized by the inelastic scattering. The inelasticity decreases with the impact 
parameter i.e. toward the peripheral region. The total differential cross-section 
is calculated according to Eq.(2.15) a,t x — 0.5, and incident energy lAGeV. The 
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prediction of the multiple scattering model, Fig. (4) shows a large peak at small values 
of g^, followed by several humps of successively decreasing height. The dashed line 
represents the behavior of the point particle scattering. It seems to tangents the 
humps of the multiple scattering model. It is expected that the two curves coincide 
at low energies where the electron wavelength is long enough and cannot probe the 
microscopic structure of the proton. Now replacing the factor that appears in 
Eq.(l.l) by the corresponding multiple scattering form as in Eq.(2.15), and 

then recalculate the enhanced structure function F,^^^. Fig. (5-7) show the model 
prediction at x = 0.08,0.125 and 0.55 respectively, compared with the prediction 
of the parton model F2 [13,14]. Appreciable deviation is observed at small values 
of X, which correspond to most deep inelastic scattering. As x increases, the two 
models come closer, where the reaction is directed toward the peripheral region that 
is characterized by single collisions only. In spite of the observed variation in the 
structure function with q^, but it still keeps its own characteristics. In other words, it 
leads to the half spin value of quarks forming the protons and also it doesn't conflict 
with the values of their fractional charge. 



f igl . eps 



Figure 1: The first three terms of the profile function for the e-p multiple scattering. The 
impulse approximation, binary collision and the third order term. 



4. Concluding Remarks. 



• The multiple scattering of the electron with constituent of the proton is very 
important to estimate its structure function. 
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FIG2.EPS 



Figure 2: The total profile funetion of the Coulomb field due to the interaction of the 
incident electron with the 3-valance quarks forming the target proton. 

• The number of scattering centers inside the proton bag plays an important role 
to explain the linear increase of the total cross-section of the electron-proton 
scattering with the center of mass energy. 

• Terms up to the double scattering shows significant shadowing effect during 
the collision. 

• The square of the scattering amplitude shows a peak followed by humps of 
decreasing height. The asymptotic behavior of which has the form of l/q4. 

• Appreciable deviation to the structure function is due to the multiple scattering 
effect, particularly for the low values of the Bjorken scaling variable x which 
correspond to the deep inelastic scattering. 
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FIG3.EPS 



Figure 3: The ratio of imaginary to real part of the total profile function of the Coulomb 
field due to the interaction of the incident electron with the 3-valance quarks forming the 
target proton. 
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FIG4 . EPS 



Figure 4: The differential cross section for the scattered electron in e-p experiment at 
14 GeV. The black line represents the prediction of the multiple scattering model and the 
dotted red one is the 1/g^ law. 



FIG5.EPS 



Figure 5: The proton structure function as deduced from the experiments SLAC-E-049 
and SLAC-E-089 at x = 0.08. The dashed line represents the modified function due to 
multiple scattering. 



9 



JHEP 2.0 DRAFT (draft, published, notoc) 



JHEP, February 1, 2008, 12:49 



FIG6.EPS 



Figure 6: The proton structure function as deduced from the experiments SLAC-E-049 
and SLAC-E-089 at x = 0.125. The dashed hne represents the modified function due to 
multiple scattering. 



FIG7.EPS 



Figure 7: The proton structure function as deduced from the experiments SLAC-E-049 
and SLAC-E-089 at x = 0.550. The dashed line represents the modihed function due to 
multiple scattering. 
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